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1.0  BACKGROUND 

This  program  was  conceived  and  has  been  carried  out  to  provide 
fundamental  insight  into  the  cascade  stalling  process  and  related  post¬ 
stall  behavior  in  the  axial-flow  compressors  of  aircraft  gas  turbine 
engines. 

The  multi-stage  axial-flow  compressor  has  undergone  continuous  and 
successful  development  by  the  aircraft  engine  community  over  the  past  40 
years.  High  pressure  ratios  with  high  efficiency  have  been  obtained, 
and  means  to  produce  adequate  stability  margin  have  been  developed  and 
employed.  Dual-spool  designs,  variable  stagger  stators,  compressor 
bleed  ports,  and  rotor  casing  tip  treatment  have  been  employed 
separately  or  in  combinations  to  increase  compressor  stability  margin 
over  the  engine  speed  range.  While  the  compressor  stability  margin, 
shown  graphically  in  Fig.  1,  has  been  the  object  of  much  interest  and 
research,  less  is  known  about  the  post-stall  operating  line. 

If  the  compressor  stability  margin  is  exceeded  as  a  result  of 
unsteady  or  off-design  engine  operation,  a  surge  or  a  rotating  stall  is 
the  result,  depending  on  the  presently  well-known  system  and  compressor 
operating  parameters.  The  compressor  can  operate  continuously  on  the 
post-stall  operating  line  if  rotating  stall  is  the  result  of  exceeding 
the  stability  margin. 
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The  rotating  stall  phenomenon  has  been  extensively  investigated 
both  analytically  and  experimentally.  The  fluid  mechanics  involved  are 
very  complex,  and  it  has  never  been  possible  to  accurately  predict  the 
exact  inception  point,  the  number  and  speed  of  the  cells,  or  the  post¬ 
stall  operating  and  recovery  line.  Improved  understanding  of  rotating 
stall  and  compressor  post-stall  behavior  through  fundamental 
experimental  investigation  was  the  goal  of  the  present  research. 

2.0  THE  RESEARCH  PROGRAM 

Prior  to  the  present  research  program,  the  Principal  Investigator 
had  carried  out  research  on  rotating  stall  and  recovery  characteristics 
employing  a  mathematical  model*.  The  model  gave  indications  that  the 
fluid  dynamic  losses  of  the  compressor  blade  row  were  of  principal 
importance  in  determining  the  post-stall  and  recovery  characteristic. 
Extensive  searches  of  the  literature  showed  that  there  was  almost  no 
data  on  the  development  of  losses  in  compressor  cascades  operated  at  the 
high  angles  of  attack  associated  with  rotating  stall.  Virtually  all 
data  sets  ended  at  an  angle  of  attack  which  produced  twice  the  minimum 
losses,  far  less  than  was  necessary  to  investigate  the  role  of  cascade 
losses  in  compressor  post-stall  behavior.  It  was  desirable  to  conduct 
additional  research  to  obtain  the  loss  data  necessary  to  evaluate  the 
predictions  of  the  model,  and  to  understand  the  developing  and  limiting 
stalling  process  in  the  compressor  cascade. 

The  program  objectives  included  the  further  development  of  the 
compressor  post-stall  mathematical  model  as  a  tool  to  understand 
the  role  of  cascade  losses  in  post-stall  behavior,  the  modification 
of  a  cascade  wind  tunnel  for  experimental  research,  and  the 
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measurement  of  cascade  losses  over  a  wide  range  of  angle  of  attack  for 
three  stagger  angles  typical  of  modern  compressor  design.  With  the  data 
and  the  improved  model,  it  was  anticipated  that  some  understanding  of 
the  underlying  fluid  mechanics  of  post-stall  behavior  could  be  deduced, 
and  that  the  hypothesis  of  the  principal  role  of  cascade  losses  in 
determining  that  behavior  could  be  checked. 

To  pursue  program  objectives,  planned  activities  were  organized 
under  three  tasks,  as  follows: 

Task  1,  Cascade  Tests  -  The  VPI  high-angle-of -attack  cascade  test 
facility  shown  in  Fig.  2  was  to  be  employed  for  experimental  investiga¬ 
tions.  Initial  tests  verified  that  accurate  cascade  loss  and  turning 
angle  data  could  be  obtained  for  flow  angles  well  in  to  the  stalling 
regime  .  Employing  the  circular  arc  airfoils  installed  in  the  facility, 
loss  and  turning  angle  data  were  to  be  obtained  for  three  stagger  angle 
settings,  25,  36.5  and  45  degrees.  The  36.5  degree  setting  corresponded 
to  the  stagger  angle  of  the  VPI  high  speed  compressor  rig,  for  which 
stalling  data  was  available.  In  addition  to  the  basic  data,  blade 
pressure  profiles  and  the  general  nature  of  the  observed  flow  separation 
were  to  be  recorded. 

Task  2,  Prediction  of  Separated  Flow  in  Cascades  -  While  the 
cascade  ..experiments  were  expected  to  produce  valuable  results,  it  was 
ultimately  desired  to  develop  prediction  techniques  for  unsteady, 
stalling  flows  in  cascades.  Starting  from  the  referenced  analysis  by  H. 
Moses  ,  it  was  planned  to  develop  a  cascade  flow  analysis  which 
accurately  predicted  post-stall  flow  behavior. 
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Fig.  2  Cascade  Geometry 
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Task  3,  Improvements  In  the  Model  for  Post-Stall  Response  In 

Coapressors  -  The  compressor  model  was  needed  to  evaluate  the 
information  gained  from  the  cascade  flow  analysis  and  experiments  for 
relevance  to  compressor  stage  performance.  Certain  improvements  to  the 
present  model  were  desirable.  These  included  improvements  in  the  method 
for  applying  the  dynamic  response  function  to  cascade  data, 
incorporation  of  a  means  for  calculating  circumferential  static  pressure 
variation,  and  the  addition  of  a  global  calculation  for  prediction  of 
compressor  flow  and  pressure  rise. 

3.0  SUMMARIZED  RESULTS 

Results  obtained  in  the  program  are  organized  according  to  the 
Tasks  listed  above. 

3.1  Cascade  Tests 

The  high-angle-of-attack  cascade  tunnel  was  suitably  configured, 
and  fundamental  experiments  to  determine  the  nature  of  the  stalling 
process  and  loss  development  were  conducted.  It  was  planned  to  measure 
the  fluid  losses  with  pneumatic  probes,  following  the  determination  of 
local  flow  directions  with  hot  film  probes.  An  extensive  program  of 
flow  visualization  was  undertaken,  so  that  the  point  measurements  of  the 
probes  could  be  properly  interpreted.  The  flow  visualization  work 
provided  some  of  the  most  important  results  of  the  investigation. 
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3.1.1  Surface  Flow  Visualization 


Oil-film  surface  flow  visualization  was  used  to  study  the  blade 

surface  flow  features  with  high-angle-of-attack  flows.  In  the  full 
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report  of  this  work  ,  approximately  50  sketches  and  photographs  describe 
the  results.  Tests  were  conducted  at  stagger  angles  of  25,  36.5,  and  45 
degrees.  A  summary  of  the  results  are  provided  here. 

Figures  3  and  4  show  photographs  of  typical  surface  flow  patterns 
produced  by  the  cascade  flows.  From  the  patterns,  observations 
concerning  the  boundary  layer  development  and  the  nature  of  stalling  in 
the  cascade  were  deduced.  A  map  summarizing  the  boundary  layer 
development  at  the  mid-span  section  of  the  cascade,  for  all  stagger 
angles  tested,  is  presented  in  Fig.  5.  This  figure  was  constructed  from 
detachment  and  reattachment  points  measured  directly  from  the  blade 
surfaces.  A  general  curve  was  fitted  through  these  points  to  divide  the 
different  boundary  layer  regions. 

At  negative  incidences,  the  blade  row  stagger  appears  to  have 
little  effect.  Laminar  flow  exists  on  the  suction  surface  for  the  first 
60  percent  chord,  whereupon  the  boundary  layer  detaches  completely. 
Separation  occurs  earlier  as  the  angle  of  attack  is  increased.  On  the 
pressure  surface,  aside  from  the  laminar  separation  bubble  at  the 
leading  edge,  the  boundary  layer  is  turbulent.  The  bubble  region 
decreases  with  increasing  angle  of  attack,  and  disappears  as  the 
incidence  becomes  positive,  at  which  the  flow  is  fully  attached.  The 
nature  of  the  boundary  layer  could  not  be  deduced  from  the  visuali¬ 
zations,  but  the  literature  suggests  it  is  mostly  laminar. 

For  small  positive  incidences,  a  laminar  separation  bubble  forms  at 


Leading  Edge 
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Leading  Edge  PRESSURE  SURFACE 


Figure  5.  Summary  of  the  Midspan  Boundary  Layer  Development  for 
Rec  =  2.0  x  105  and  Tu  =  1.1  percent 


the  leading  edge  of  the  suction  surface  and  remains  present  until  stall 
is  complete.  The  reattachment  data  show  much  scatter.  It  is  felt  that 
the  bubble  length  did  not  actually  change  much  with  respect  to  angle  of 
attack  or  stagger  angle,  and  that  the  scatter  is  due  to  an  error  in 
measurement.  The  reattachment  line  in  most  cases  was  difficult  to 
distinguish.  Moreover,  it  is  believed  that  the  appearance  of  the 
leading  edge  separation  bubble  on  the  pressure  surface  at  negative 
incidences  and  on  the  suction  surface  at  positive  incidences  is  due  to 
the  relative  sharpness  of  the  blade  nose,  and  that  its  effect  on  the 
overall  features  of  stall  is  insignificant. 

The  change  in  the  blade  surface  flow  characteristics  as  the  angle 
of  attack  is  varied  from  negative  to  positive  incidence  is  believed  to 
occur  suddenly.  No  physical  evidence  was  obtained  in  this  work  to  fix 
the  incidence  of  transition.  However,  as  shown  in  the  figure,  it  is 
reasonable  to  assume  that  transition  occurs  at  zero  incidence. 

The  predominant  characteristic  of  the  suction  surface  stall  at 
positive  incidence  is  the  rapid  progress  of  the  turbulent  separation 
zone  toward  the  leading  edge  of  the  airfoil.  The  effect  of  the  stagger 
angle  is  manifested  in  the  severity  of  stall  at  equal  angles  of 
attack.  The  higher  stagger  cascade  has  a  narrower  operating  range.  At 
midspan,-  full  stall  occurred  at  an  angle  of  attack  near  20  degrees  for 
the  45  degree  stagger  setting,  and  near  22.5  degrees  for  the  two  other 
staggers. 

The  general  nature  of  the  observed  three-dimensional  stalling 
process  and  surface  flow  under  fully  stalled  conditions  is  summarized  in 
schematic  form  in  Fig.  6.  Preceding  and  following  the  inception  of 


Angle  of  Attack, 


complete  separation  at  the  leading  edge,  the  flow  at  the  blade  suction 
surface  is  essentially  two-dimensional.  Corner  stall  is  the  predominant 
mechanism  for  three-dimensionality  preceding  full  stall.  Corner  stall 
appears  at  small  positive  incidence  and  becomes  physically  more  severe 
with  increases  in  the  angle  of  attack.  At  operation  near  full  stall, 
the  turbulent  separation  was  delayed  by  the  corner  stall,  but  because  of 
the  high  aspect  ratio  blades  used,  this  effect  was  concentrated  in 
regions  away  from  midspan.  It  is  easy  to  imagine  from  these  results  how 
the  turbulent  separation  in  the  midspan  region  would  have  been  altered 
and  delayed  if  a  lower  aspect  ratio  cascade  had  been  employed  in  this 
investigation. 

The  observed  boundary  layer  behavior  was  related  to  axial  velocity 
ratio  (AVR)  measurements.  It  was  found  that  an  increase  in  the  angle  of 
attack  resulted  in  a  rise  in  the  AVR.  The  rate  of  increase  slowed  in 
the  fully  stalled  flow  regime.  The  AVR  data  showed  a  maximum  value  of 
about  1.15.  The  rise  in  the  AVR  above  unity  paralleled  the  growth  seen 
in  the  corner  stall  region.  In  the  fully  stalled  flow  regime,  the 
diffusion  by  the  cascade  levels,  as  does  the  pressure  rise,  thereby 
limiting  any  further  increase  in  the  corner  stall  region  and  the  AVR. 
Here,  the  surface  streamlines  were  found  to  be  straight  for  a  large 
percentage  of  the  blade  height,  but  the  AVR  did  not  return  to  unity. 
This  illustrates  the  possible  error  in  describing  tests  as  being  truly 
two-dimensional  based  solely  on  surface  flow  visualization  results. 
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3.1.2  Reynolds  Number  Effects 

Additional  blade  surface  flow  visualization  tests  were  conducted  to 
provide  physical  evidence  as  to  whether  the  test  Reynolds  number  of  2.0 
x  10^  was  above  the  critical  Reynolds  number  range.  This  series  of 
experiments  was  carried  out  with  the  blade  cascade  set  at  a  stagger  of 
36.5  degrees.  Visualizations  were  obtained  for  Reynolds  number  of  2.0  x 
105  ,  1.5  x  10^,  and  1.0  x  10^  at  angles  of  attack  of  5,  15,  20,  25,  and 
35  degrees.  The  Reynolds  number  was  varied  by  adjusting  the  inlet 
velocity. 

In  summary,  the  rapid  deterioration  in  the  performance  of  a 
compressor  cascade  with  decreasing  Reynolds  number  is  generally  thought 
to  result  from  the  appearance  and  growth  of  laminar  separation  bubbles 
on  the  suction  surface  of  the  blades.  The  Reynolds  number  test  results 
showed  relatively  small  changes  in  boundary  layer  development  with 
decreasing  Reynolds  number,  and  any  increase  in  losses  is  believed  to  be 
minimal.  The  pressure  surface  flow  showed  no  variation  with  Reynolds 
number.  On  the  suction  surface,  some  small  increases  occurred  in  the 
size  of  the  laminar  and  turbulent  separation  zones.  The  leading  edge 
separation  bubble  observed  on  the  pressure  surface  at  negative 
incidences,  and  on  the  suction  surface  at  positive  incidences,  is 
believed  to  be  due  to  the  relative  sharpness  of  the  blade  nose  and  is 
not  a  Reynolds  number  effect.  The  only  real  doubt  as  to  whether  the 
Reynolds  number  of  2.0  x  10^  is  above  the  critical  range  involves 
laminar  separation  observed  on  the  suction  surface  at  negative 
incidences.  It  is  possible  that  at  a  higher  Reynolds  number,  only  a 
partial  separation  would  be  obtained  with  turbulent  reattachment. 
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3.1.3  Tip  Clearance  Effects 


Surface  flow  visualization  was  used  to  investigate  the  effects  of 
tip  clearance  on  boundary  layer  development  across  the  blade  span.  For 
this  investigation,  bridge  pieces  were  removed  from  the  test  section 
leaving  a  clearance  of  gap/chord  ratio  equal  to  4.8  percent  between  the 
left  sidewall  and  blade  tips.  With  the  blade  cascade  set  at  a  stagger 
angle  of  36.5  degrees,  tests  were  performed  at  angles  of  attack  of  15, 
20,  25,  and  30  degrees.  Sample  sketches  of  the  blade  surface  flow 
patterns  obtained  with  and  without  tip  clearance  are  shown  in  Fig.  7. 
In  the  figure,  the  gap  end  is  to  the  left.  Additional  photos  and 
sketches  are  shown  in  Ref.  3. 

In  summary,  at  positive  incidences  prior  to  full  stall,  tip  leakage 
was  very  effective  in  preventing  the  premature  separation  on  the  blade 
surface  associated  with  corner  stall.  Separation  was  actually  delayed 
for  a  short  distance  from  the  gap.  This  distance  increased  somewhat 
with  increasing  angle  of  attack.  Inward  of  the  region  under  the  direct 
influence  of  the  leakage  flow,  the  limiting  streamlines  converged 
towards  midspan,  but  at  an  angle  somewhat  less  than  that  caused  by 
corner  stall  with  zero  clearance.  For  operation  in  the  fully-stalled 
region,  the  gap  produced  some  unstalled  flow  near  the  gap,  but  for  all 
practical  purposes  the  effect  was  considered  small.  Further,  a  three- 
dimensional  flow  resulted  at  the  gap  end  which  was  equal  in  extent  to 
that  found  at  the  zero  clearance  end. 


Separation  Bubble  Leading  Edge 

Sucuon  Side 


(a)  Zero  Clearance 


Leakage  Induced 


(b)  Clcarance/Chord  =  4.8  percent 


Figure  7.  Sketches  of  the  Surface  Flow  Patterns  Showing  Tip  Clearance  Effects 
for  £  =  36.5°  and  a*  =  15° 


3.1.4  Conclusions,  Surface  Flow  Visualization 

Oil-flow  visualizations  were  conducted  in  a  solid  wall,  high  aspect 
ratio,  compressor  cascade  for  the  purpose  of  describing  the  transition 
from  corner  stall  to  full  blade  stall,  and  the  blade  surface  flow  under 
fully  stalled  conditions.  Results  were  obtained  for  stagger  angles  of 
25,  36.5,  and  45  degrees. 

Similar  trends  in  boundary  layer  development  were  exhibited  for  the 
three  stagger  angles.  At  negative  incidences,  the  suction  surface 
showed  substantial  areas  of  laminar  flow.  The  laminar  boundary  layer 
terminated  in  a  complete  separation  of  flow  from  the  blade  surface.  On 
the  pressure  surface  a  laminar  separation  bubble  existed  at  the  leading 
edge  with  the  remaining  boundary  layer  turbulent. 

At  small  positive  incidences  a  laminar  separation  bubble  developed 
at  the  leading  edge  of  the  suction  surface,  followed  by  turbulent 
reattachment  a  short  distance  downstream.  The  turbulent  boundary  layer 
separated  near  the  trailing  edge.  As  the  angle  of  attack  was  increased 
to  the  point  of  a  complete  leading  edge  separation,  the  turbulent 
separation  progressed  upstream,  with  no  large  changes  in  the  laminar 
bubble.  The  boundary  layer  on  the  pressure  surface  was  fully  attached 
for  the  range  of  positive  incidences  tested;  the  boundary  layer 
structure  could  not  be  determined  from  the  experiments. 

The  effects  of  stagger  were  revealed  in  the  severity  of  the  stall 
at  positive  incidences.  Transition  to  complete  separation  at  the 
leading  edge  occurred  at  lower  angles  of  attack  for  the  higher  stagger 
angles. 

The  most  interesting  discoveries  from  the  work  deal  with  the 


profile  boundary  layer  three-dimensionality.  Geometrically  speaking, 
the  flow  on  the  pressure  surface  was  found  to  be  virtually  two- 
dimensional  for  all  tested  angles  of  attack.  The  suction  surface  flow 
was  shown  to  be  essentially  two-dimensional  preceding  and  following  the 
transition  to  a  fully  separated  flow  at  the  leading  edge.  Corner  stall 
was  observed  at  small  positive  incidence  and  grew  in  severity  with 
increasing  angles  of  attack.  Near  full  stall  at  midspan,  turbulent 
separation  was  delayed  in  regions  adjacent  to  the  corner  stall  zones. 
For  fully-stalled  conditions,  the  three-dimensional  mechanism  took  the 
form  of  recirculating  flow  regions  at  the  ends  of  the  blades. 

Special  oil-flow  tests  were  conducted  to  study  Reynolds  number  and 
tip  clearance  effects  on  the  blade  boundary  layer  development  under 
normal  and  fully-stalled  conditions.  No  pronounced  changes  in  boundary 
layer  development  were  observed  in  the  Reynolds  number  range  from  1.0  x 
105  to  2.0  x  10^.  The  tip  clearance  tests  showed  for  this  cascade  that 
the  gap/chord  ratio  of  4.8  prevented  corner  stall  and  actually  delayed 
turbulent  separation  for  some  distance  from  the  gap.  After  the  flow  was 
fully  separated,  the  gap  had  only  small  effects  concentrated  near  the 
end  of  the  blade. 

3.1.5  Smoke  Flow  Visualization  Results 

Smoke  flow  visualization  tests  were  conducted  to  qualitatively 
examine  the  nature  of  the  stalled  flow  and  to  determine  the  existence 
and  type  of  unsteadiness.  Separated  flows  are  usually  unsteady. 
However,  an  ordered  unsteadiness  in  the  form  of  a  propagating  stall 
could  also  exist  when  testing  cascades  of  airfoils.  The  existence  and 


type  of  unsteadiness  is  important  in  itself,  particularly  if  design 
variables  affect  the  occurrence  of  propagating  stall  and  the  propagating 
speed.  It  was  also  important  to  know  if  the  flow  was  unsteady  when 
evaluating  the  data  from  the  pressure  and  velocity  measurements. 

Still  shot  photographs  were  intially  taken  over  a  range  of  angle  of 
attack  for  the  cascade  with  a  36.5  degree  stagger.  These  photographs 
were  taken  to  determine  if  high  speed  photography  would  reveal  the  flow 
structures  of  interest.  With  the  high  velocity  of  the  flow,  the  unaided 
eye  could  not,  at  the  very  least,  resolve  the  separated  regions  and 
obviously  could  not  identify  the  unsteadiness.  The  separated  flow 
regions  were  visible  when  strobe  lights  were  used  to  illuminate  and  slow 
the  flow.  However,  the  unsteadiness  still  could  not  be  observed. 
Figures  8,  9,  and  10  are  photographs  of  the  flow  at  15,  25,  and  35 
degrees  angles  of  attack.  The  flow  in  the  photographs  was  frozen  by 
using  seven  Strobotacs  to  provide  the  flash.  These  strobes  lights  were 
connected  in  parallel  and  triggered  by  the  camera. 

In  Fig.  8  the  smoke  reveals  the  free  stream  turbulence  and  a 
possible  small  separated  region  near  the  trailing  edge  of  the  airfoil. 
For  this  angle  of  attack  the  boundary  layers  appear  to  be  attached  over 
most  of  the  blade.  Mixing  in  the  wake  of  the  airfoils  is  evident  from 
the  more,  diffuse  smoke  in  this  region. 

Figures  9  and  10  show  the  flow  for  angles  of  attack  of  25  and  35 
degrees,  respectively.  For  both  of  these  photographs,  large  separated 
regions  are  evident.  It  appears  that  the  flow  separates  from  the 
leading  edge  of  the  blades  and  that  the  size  of  the  separated  region 
increases  (as  expected)  with  an  increase  in  angle  of  attack.  It  is  also 


Figure  10.  Photograph  of  the  smoke  flow  visualization  for  C  =  36.5°  and  a  =  35 
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evident  in  both  Figs.  9  and  10  that  the  separated  regions  at  this 
instant  in  time  are  not  identical  in  the  two  neighboring  passages.  In 
Fig.  9  the  separated  region  in  the  middle  passage  is  larger  than  the 
separated  region  in  the  passage  below,  and  has  diffuse  smoke  in  the 
entire  region.  In  the  passage  below  the  middle  passage,  the  separated 
region  is  smaller.  The  differences  observed  in  the  two  passages  of  Fig. 
9  are  similar  to  what  would  be  expected  if  propagating  stall  was 
present.  However,  conclusive  evidence  with  regard  to  the  existence  of 
propagating  stall  would  require  sequential  photographs  taken  at  very 
short  intervals  (or  high  speed  motion  pictures). 

The  still  shot  photographs  verified  that  the  structure  of  the  flow 
could  be  observed  with  the  present  smoke  visualization  system.  However, 
it  also  became  apparent  that  high  speed  movies  would  be  necessary  to 
determine  if  propagating  stall  occurs.  A  Hycam  II  high  speed  rotating 
prism  camera  was  used  to  photograph  the  flow  at  4000  frames  per  second. 
Motion  pictures  of  the  flow  were  created  for  test  cases  covering  a  full 
range  of  angle  of  attack  and  for  staggers  of  25,  36.5,  and  45  degrees. 
For  selected  cases,  film  was  taken  for  high  and  low  inle.  velocities 
corresponding  to  Reynolds  numbers  of  approximately  200,000  and  70,000. 

Viewing  the  films  at  a  normal  projection  speed  of  16  frames/second 
slowed  the  flow  by  a  factor  of  250  to  1.  For  the  25  degree  stagger 
cascade,  propagating  stall  was  not  observed  even  at  angles  of  attack  as 
large  as  45  degrees.  For  the  cascades  with  staggers  of  36.5  and  45 
degrees,  propagating  stall  was  observed  at  all  angles  of  attack  of  20 
degrees  and  greater.  Contrary  to  what  may  have  been  expected,  the  flow 
did  not  restabilize  at  high  angles  of  attack,  but  continued  to  reveal 
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propagating  stall  superimposed  on  a  mean  full  stall  from  the  leading 
edge.  The  propagating  stall  evident  in  the  films  showed  the  classical 
blockage  of  a  passage,  the  flow  diversion  upstream  of  the  cascade  around 
the  blockage,  and  the  subsequent  recovery  (or  complete  flushing  out)  in 
the  passage  behind  the  propagating  cell.  Ahead  of  the  stall  cell,  the 
increase  in  angle  of  attack  caused  by  the  flow  diverting  around  the 
blockage  resulted  in  the  next  blade  stalling.  It  was  also  observed  that 
stall  propagation  was  associated  with  almost  complete  blockage  of  a 
passage  and  not  just  an  increase  in  blockage  from  the  mean.  In  the  25 
degree  stagger  cascade,  the  stall  did  not  propagate,  and  in  this  case 
the  passage  was  observed  to  never  be  completely  blocked. 

Piatt^*  examined  the  films  which  were  created  in  this  study  in 
detail.  He  determined  the  number  of  occurrences  of  stall  propagation 
and  the  propagation  speeds  for  each  test  case  with  the  36.5  and  45 
degree  stagger  cascade.  The  complete  results  of  his  study  are  presented 
in  Ref.  4,  along  with  groups  of  sequential  frames  from  the  films  to 
provide  supporting  evidence.  Figure  11  is  a  summary  of  Piatt's  obser¬ 
vations  of  the  history  of  a  typical  propagating  stall  event.  He 
concluded  that  the  propagation  speed  was  nearly  proportional  to  the 
inlet  velocity.  A  small  decrease  in  the  propagation  speed  was  observed 
at  20  degrees  of  angle  of  attack,  which  was  associated  with  the 
occurrence  of  slightly  less  than  complete  blockage  of  the  passage.  Once 
an  angle  of  attack  was  reached  where  the  passage  became  fully  blocked, 
no  further  changes  in  the  propagation  speed  occurred. 
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re  11.  Observed  Sequence  of  Events  in  Propagating  Stall 


3.1.6  Test  Section  Inlet  Mean  Velocity  Measurements 


Russ  [3]  conducted  surveys  of  the  velocity  field  upstream  of  the 
cascade  for  the  purpose  of  evaluating  the  uniformity  of  the  flow. 
Surveys  were  made  from  side  wall  to  side  wall  in  the  spanwise  direction 
and  over  a  distance  of  4  blade  spacings  in  the  pitch  direction.  Russ' 
measurements  showed  that  the  flow  in  the  spanwise  direction  is  uniform 
within  1  m/s  outside  of  the  wall  boundary  layers.  The  boundary  layers 
on  the  side  walls  vary  in  thickness  from  approximately  14  mm  to  29  mm, 
depending  on  which  wall  is  considered  and  the  pitch  location  of  the 
measurement.  Even  for  the  maximum  measured  boundary  layer  thickness  of 
29  mm,  the  blades,  which  have  a  span  of  308.0  mm,  are  subjected  to 
uniform  flow  over  81%  of  the  span.  Thus,  highly  uniform  inlet  flow  in 
the  span  direction  was  present  during  the  tests. 

Russ  conducted  upstream  pitchwise  traverses  at  three  span  locations 
and  three  Reynolds  numbers.  All  traverses  exhibited  an  approximately 
sinusoidal  variation  in  the  velocity,  with  a  wavelength  equal  to  the 
blade  spacing.  Since  the  traverses  were  conducted  at  a  25  degree  angle 
of  attack,  the  upstream  variation  in  velocity  can  be  attributed  to  the 
divergence  around  the  regions  blocked  by  stalled  flow.  The  peak 
amplitude  of  the  velocity  variation  had  its  greatest  value  of 
approximately  3%  of  the  mean  velocity  for  the  tests  with  the  largest 
inlet  velocity. 

Another  means  of  evaluating  the  periodicity  of  the  flow  is  by 
measuring  the  pressure  from  the  static  wall  taps  located  upstream  of  the 
cascade.  These  taps  were  spaced  at  intervals  equal  to  one  blade 
spacing,  and  thus  periodic  flow  would  be  indicated  by  equal  pressures. 
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Measurements  of  the  pressures  from  these  taps  were  made  during  the 
pressure  measurement  series  of  tests.  A  variation  in  static  pressure 
from  the  top  to  the  bottom  does  exist  for  every  nonzero  cascade  inlet 
angle,  a^.  The  variation  is  simply  the  result  of  the  lower  taps  being 
at  a  greater  distance  downstream  of  the  nozzle,  where  a  larger  pressure 


drop  has  occurred  due  to  viscous  friction  in  the  boundary  layers. 
Because  of  the  source  of  the  variation,  it  would  be  expected  to  have  the 
greatest  variation  at  the  higher  inlet  angles.  For  the  maximum  inlet 
angle  of  a ^  =  75  degrees,  the  difference  in  static  pressure  between  the 
taps  located  five  spacings  apart  is  11.2  mm  l^O.  This  corresponded  to 
only  8.5%  of  the  inlet  dynamic  head  for  this  worst  case.  At  the  minimum 
inlet  angle  tested,  a j  =  25  degrees,  the  difference  in  static  pressure 
over  three  blade  spacings  was  3.49  mm  1^0,  which  was  3%  of  the  inlet 
dynamic  head.  For  the  purpose  of  modeling  two-dimensional  cascade  flow, 
the  inlet  velocity  field  was  therefore  shown  to  be  uniform  in  the 
spanwise  direction  and  essentially  periodic  in  the  pitchwise  direction. 

3.1.7  Turbulence  Intensity  Measurements 

Free  stream  turbulence  intensity  levels  are  known  to  affect  the 
boundary  layer  transition,  the  type  and  location  of  separation  and, 
consequently,  the  performance  of  a  compressor  blade  row.  Although  it  is 
not  necessary  to  minimize  the  freestream  turbulence,  it  is  important  to 


measure  the  turbulence  in  order  to  completely  characterize  the  inlet 
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conditions.  The  turbulence  intensity,  — — ,  was  measured  at  an  array  of 

u 

points  for  four  difference  inlet  velocities  corresponding  to  Reynolds 
number  based  on  the  blade  chord  of  200000,  150000,  100000,  and  57000. 


The  array  of  measurement  locations  consisted  of  three  span  locations  of 


z  =  79,  156,  and  232  mm.  At  each  span  position,  measurements  were  made 
at  25.4  mm  increments  in  the  pitch  direction.  These  data  were  obtained 
with  a  cascade  stagger  angle  of  45  degrees  and  an  angle  of  attack  of  10 
degrees.  Preliminary  tests  at  an  angle  of  attack  of  25  degrees 
exhibited  severe  unsteadiness,  probably  as  a  result  of  propagating 
stall. 

The  turbulence  intensities  ranged  from  1.  I  to  1.3%  for  all  inlet 
velocities  and  measurement  locations.  The  variation  in  the  intensity 
was  relatively  small  and  was  considered  insignificant.  The  point  to 
point  variation  in  mean  velocity  was  less  than  1  m/s. 

3.1.8  Velocity  Measurements  and  Numerical  Predictions 

The  magnitude  and  angles  of  the  flow  velocity  in  the  cascade  were 
obtained  experimentally  from  measurements  with  the  dual  split  film 
probe.  The  time  mean  flow  angles  downstream  of  the  cascade  were  the 
main  quantities  of  interest,  being  required  to  align  the  pneumatic  probe 
with  the  flow  in  the  subsequent  experiments.  For  angles  of  attack  of 
15,  20,  25,  and  30  degrees,  surveys  across  the  passage  at  approximately 
20%  and  50%  chord  were  also  conducted  to  further  define  the  extent  of 
the  separated  region  within  the  blade  passage.  The  four  signals  from 
the  hot  film  probe  were  processed  to  produce  a  velocity  measurement,  and 
twenty  measurements  were  averaged  at  each  point.  Plots  of  the  averaged 
velocity  vectors  provide  a  pictorial  indication  of  what  the  flow  is  like 
through  the  cascade.  Figures  12  through  23  show  the  measured  velocity 


vectors.  All  the  data  presented  in  these  figures  is  for  a  nominal 
Reynolds  number  of  200000. 

For  a  given  stagger  angle,  the  variation  of  the  flow  fields  with 
angle  of  attack  are  consistent  with  the  expected  trends.  At  the  lower 
angle  of  attack  of  15  degrees,  Figs.  12,  16,  and  20  show  no  reverse  flow 
indicated  by  the  measured  vectors,  and  the  vectors  increase  in  magnitude 
as  the  suction  surface  is  approached.  The  acceleration  of  the  flow  near 
the  suction  surface  is  a  potential  flow  effect  which  results  in  the 
creation  of  lift  on  the  airfoil.  As  the  angle  of  attack  is  increased  to 
20  degrees,  Figs.  12,  17,  and  21  show  that  viscous  effects  near  the 

suction  surface  have  increased,  resulting  in  a  much  larger  boundary 
layer  on  the  blades  and  a  larger  wake  downstream.  At  this  angle  of 
attack,  no  reverse  flow  is  indicated  by  the  measurements,  although 
separated  flow  was  indicated  to  some  extent  by  the  surface  flow 
visualization.  Therefore,  it  can  be  assumed  that  measurements  at  this 
angle  of  attack  could  not  be  made  near  enough  to  the  surface  to  resolve 
the  separated  flow  region. 

The  results  for  angles  of  attack  of  25  and  30  degrees  are  given  in 

Figs.  14  and  15,  18  and  19,  and  22  and  23.  At  these  angles  of  attack 

the  measurements  for  all  three  stagger  angles  indicate  reversed  flow. 
The  separated  flow  regions  are  seen  to  increase  in  size  as  the  angle  of 
attack  is  increased  from  25  to  30  degrees. 

The  most  significant  observation  to  be  made  from  Figs.  12  through 
23  is  the  change  in  the  nature  of  the  flow  associated  with  the  changes 
in  stagger  angle.  For  the  cascade  with  a  stagger  of  25  degrees,  Figs. 

14  and  15  show  that  two  distinct  regions  exist  within  the  blade 
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Figure  18.  Me 


Figure  19.  Measured  Velocities  in  Blade  Passage  Region 
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Figure  21.  Measured  Velocities  in  Blade  Passage  Reg 
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Figure  23.  Measured  Velocities  in  Blade  Passage  Region 
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have  the  same  general  shapes  and  exhibit  the  same  trends  with  respect  to 
stagger  as  do  the  measured  velocity  profiles.  Figure  24  shows  a  sample 
comparison  of  the  numerically-predicted  velocity  field  with  the  measured 
values,  for  an  angle  of  attack  of  30  degrees.  The  obvious  reasons  for 
the  differences  in  the  absolute  values  of  the  velocities  are  the  approx¬ 
imations  and  simplifying  assumptions  made  in  the  numerical  prediction. 
Even  at  an  angle  of  attack  of  15  degrees,  the  numerical  results  indicate 
some  reverse  flow.  This  separation  is  most  likely  the  result  of  the 
sharp  leading  edge  associated  with  approximating  the  blades  as  flat 
plates.  The  extent  of  the  reverse  flow  is  predicted  by  the  analysis  to 
be  larger  for  the  higher  stagger  angles.  This  result  is  consistent  with 
the  observation  from  the  surface  flow  visualization  that  full  stall 
occurs  at  a  slightly  lower  angle  of  attack  for  cascades  with  greater 
stagger. 

For  the  larger  angles  of  attack,  the  velocity  of  the  reverse  flow 
predicted  numerically  was  found  to  have  greater  magnitude,  but  with  the 
region  of  the  separated  flow  extending  a  smaller  distance  into  the 
passage  than  the  corresponding  measured  flow.  This  difference  is 
thought  to  be  due  primarily  to  the  curvature  of  the  suction  surface  of 
the  real  blade  rather  than  due  to  the  finite  thickness  of  the  real 
blade.  Full  details  of  the  velocity  measurements  and  theory  development 
are  provided  in  Ref.  5. 

3.1.9  Cascade  Pressure  Measurements 

Figure  25  shows  an  example  of  total  pressure  and  static  pressure 
losses  measured  behind  the  25  degree  stagger  cascade,  for  a  very  high 
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angle  of  attack  of  45  degrees.  Predictions  of  the  previously-mentioned 


theory  are  shown  for  comparison.  At  this  angle  of  attack,  total 
pressure  losses  are  extremely  high,  and  extend  over  a  large  percentage 
of  the  cascade  pitch.  Large  static  pressure  variations  are  measured. 
Agreement  between  theory  and  experiment  is  good  for  the  lower  loss 
regions,  but  the  theory  significantly  under-predicts  the  high  levels  of 
static  and  total  pressure  coefficient.  Much  better  agreement  is 
obtained  for  lower  angles  of  attack.  Complete  measurements  and 
comparisons  are  presented  in  Ref.  5. 

Pressure  measurements  were  made  employing  static  taps  on  the 
suction  and  pressure  surfaces  of  the  cascade  blades.  The  taps  were 
located  at  mid-span,  and  in  the  centermost  blade  passage.  Measured  and 
predicted  pressure  coefficients  are  shown  in  Fig.  26,  for  the  same  flow 
conditions  shown  in  Fig.  25.  Here,  agreement  is  excellent,  especially 
considering  the  extremely  high  angle  of  attack  of  the  flow. 

3. 1. 10  Cascade  Performance 

The  ultimate  goal  of  the  cascade  measurement  program  was  to  provide 
improved  understanding  of  the  cascade  stalling  process,  and  to  provide 
measurements  of  the  fluid  dynamic  losses  of  the  cascade  at  high  angle  of 
attack  for  three  stagger  angles.  The  previously-described  measurements 
were  utilized  to  compute  the  cascade  losses  with  results  shown  in  Fig. 
27.  The  loss  parameter  is  the  mass-averaged  total  pressure  loss 


coefficient  immediately  behind  the  cascade,  referenced  to  the  mean 
velocity  entering  the  cascade. 

A  principal  feature  of  the  loss  behavior  is  the  lower  rate-of-rise 
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Figure  25.  Measured  and  Predicted  Cascade  Losses 
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Figure  27.  Flow  Losses  for  Three  Stagger  Angles 
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of  the  losses  for  25  degrees  stagger,  as  compared  with  that  for  the  36.5 
and  45  degree  stagger  cases.  This  result  is  assumed  to  be  due  to  the 
observed  differences  in  cascade  blockage  for  the  three  stagger  angles, 
as  previously  discussed.  The  results  are  very  important  with  regard  to 
the  model  for  rotating  stall,  which  has  been  found  to  be  primarily 
influenced  by  the  slope  of  the  loss  curve  after  stalling  has  occurred. 

The  blade  normal  force  coefficient  was  calculated  from  the  measured 
velocity  and  pressure  profiles,  with  the  result  shown  in  Fig.  28. 
Again,  a  significant  difference  is  noted  for  the  25  degree  stagger  angle 
results. 

Loss  and  normal  force  coefficient  predictions  of  the  mathematical 
model  are  shown  in  Figs.  29  and  30.  Trends  are  correctly  predicted,  but 
differences  between  36.5  and  45  degree  stagger  angle  predictions  are 
larger  than  experimentally  observed.  Reasons  for  the  differences  are 
being  assessed. 

3.1.11  Conclusions,  Cascade  Measurements 

Extensive  measurements  of  cascade  flow  characteristics  and 
performance  were  accomplished.  Results  were  obtained  for  three  stagger 
angles,  as  planned.  Significant  effects  of  stagger  angle  were  noted, 
which  are  expected  to  influence  predictions  of  the  post-stall  compressor 
stage  model.  Supporting  flow  visualizations  provided  new  insight  into 
the  cascade  stalling  process.  The  data  obtained  were  compared  with 


predictions  of  a  mathematical  model  developed  under  the  program,  with 
good  results.  For  future  use,  all  data  obtained  are  reported  in  Refs. 
3,  4,  and  5. 
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Figure  23.  Normal  Force  Coefficient  for  Cascade  Blade  from  Measurements 
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Figure  30.  Calculated  Normal  Force  Coefficient  for  Cascade 


3.2  Prediction  of  Separated  Flow  in  Cascades 


A  mathematical  model  for  prediction  of  separated  flows  in  cascades 
was  developed  under  the  program.  As  previously  noted,  the  flow 
simulation  was  obtained  by  numerically  solving  the  two-dimensional, 
steady,  constant,  viscosity,  incompressible  Navier-Stok.es  equations 
using  the  SIMPLER  algorithm  of  Patankar.  Comparisons  of  the  theory 
predictions  and  measured  results  were  shown  in  the  previous  sections. 
The  complete  development  of  the  theory  is  given  in  Ref.  5. 

3.3  I tnr  •" o  ve me n t s  in  Post-Stall  Compressor  Stage  Model 

At  the  beginning  of  the  program,  a  post-stall  predictive  model 
existed  from  previous  work*,  but  improvements  in  the  model  were 
desirable.  A  description  of  the  model  and  progress  made  was  published^, 
and  a  copy  of  the  paper  is  included  here  as  Appendix  A.  As  a  result  of 
work  under  the  present  program,  the  model  is  presently  capable  of 
predicting  a  compressor  stage  post-stall  characteristic,  and  accomodates 
a  non-uniform  downstream  pressure  boundary  behind  the  simulated  rotor. 
Work  continues  to  add  a  stator  row  to  the  model,  thus  simulating  a  full 
stage.  With  the  model  and  the  data  obtained  from  the  investigation,  new 
insights  have  been  provided,  and  work  continues  to  add  to  the 
understanding  of  the  mechanisms  of  compressor  post-stall  behavior. 

4.0  PUBLICATIONS 

1.  Cousins,  W.  T.  ,  and  W.  F.  O'Brien,  "Compressor  Stage  Post-Stall 
Analysis,"  AIAA  Paper  No.  75-1349. 

When  the  addition  of  a  downstream  stator  row  model  is  complete,  it 
is  planned  to  submit  a  revised  version  of  this  paper  to  the  AIAA 
Journal  of  Propulsion  and  Power. 
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ABSTRACT 

This  study  was  conducted  for  the  purposes  of  providing  a  more 
fundamental  understanding  of  the  flow  and  to  provide  performance  data 
for  fully  stalled  cascades.  Performance  data  for  fully  stalled 
compressor  blade  rows  are  required  to  model  compressor  stall  phenomena 
using  currently  available  techniques.  Cascades  of  a  single  blade 
geometry  and  a  solidity  of  unity  were  studied  for  three  stagger  angles 
and  over  the  full  range  of  angle  of  attack  extending  well  into  the 
stalled  flow  regime.  The  Reynolds  number  was  also  varied  for  a  limited 
number  of  cases. 

Part  I  of  this  paper  presents  the  results  obtained  from  surface 
flow  visualization  tests,  smoke  flow  visualization  tests,  and  time-mean 


velocity 

measurements. 

Surface 

flow 

visualization 

studies 

using 

the 

oil-film 

visualization 

technique 

we  re 

conducted 

to 

examine 

the 

two- 

dimensionality  of  the 

flow  and 

to 

determine 

the 

location  of 

the 

separation  line  on  the  blades.  Smoke  flow  visualization  studies  were 
conducted  for  determining  the  existence  and  type  of  unsteadiness.  The 
visualized  flow  was  filmed  at  a  rate  of  4000  frames  per  second  using  a 
rotating  prism  high  speed  motion  picture  camera,  and  was  subsequently 
studied  in  detail  by  examining  the  films  at  a  normal  projection  speed  of 
L6  frames  per  second.  The  surface  flow  visualization  and  the  smoke  flow 
visualization  showed  that  the  blade  stagger  angle  is  a  key  parameter  in 
determining  the  location  of  the  separation  line  and  the  occurrence  of 
propagating  stall.  Time  mean  velocity  measurements  obtained  with  a  dual 
hot  split  film  probe  also  showed  that  the  separated  velocity  profiles 
within  the  blade  passages  and  the  downstream  profiles  in  the  wake  have 
distinctly  different  characteristics  depending  on  the  stagger  angle. 
The  characteristics  of  the  velocity  profiles  can  be  related  by 
physically  reasoning  to  the  occurrence  of  propagating  stall  and  the 
performance  of  the  fully  stalled  cascade,  which  is  presented  in  Part  II 
of  this  paper. 

Separated  Flow  in  a  Compressor  Cascade  -  Part  II: 

Cascade  Performance 

ABSTRACT 

This  study  was  conducted  for  the  purposes  of  providing  a  more 
fundamental  understanding  of  the  flow  and  to  provide  performance  data 
for  fully  stalled  cascades.  Performance  data  for  fully  stalled 
compressor  blade  rows  are  required  to  model  compressor  stall  phenomena 
using  currently  available  techniques.  Cascades  of  a  single  blade 
geometry  and  a  solidity  of  unity  were  studied  experimentally  for  three 


stagger  angles  and  over  the  full  range  of  angle  of  attack,  extending  well 
into  the  stall  flow  regime.  The  Reynolds  number  was  also  varied  for  a 
limited  number  of  cases. 

Part  II  of  this  paper  presents  the  results  from  the  velocity  and 
pressure  measurements  made  in  the  cascade  and  the  overall  cascade 
performance  evaluated  from  these  measurements.  In  addition,  results 
from  a  numerical  simulation  of  the  flow  through  a  cascade  of  flat  plate 
airfoils  is  used  to  further  illustrate  the  effects  of  blade  stagger  and 
to  define  the  correct  limits  for  the  cascade  performance  as  an  inlet 
angle  of  90  degrees  is  approached.  The  flow  simulation  was  obtained  by 
numerically  solving  the  two-dimensional,  steady,  constant  viscosity, 
incompressible  Navier-Stokes  equations  using  the  SIMPLER  algorithm  of 
Patankar.  The  cascade  performance  is  presented  in  terms  of  the  normal 
force  coefficient  and  the  mass  averaged  total  pressure  loss  through  the 
cascade  as  a  function  of  angle  of  attack.  The  results  indicate  that  the 
slope  of  the  total  pressure  loss  curve  for  the  flow  immediately 
downstream  of  the  cascade  is  steeper  for  cascades  with  greater  stagger, 
and  the  losses  continue  to  rise  as  the  limiting  case  of  90  degree  inlet 
flow  is  approached.  The  losses  calculated  for  the  flow  when  it  becomes 
fully  mixed  are  not  as  greatly  influenced  by  stagger.  The  normal  force 
coefficient  was  found  to  increase  to  a  peak  value  near  the  angle  of 
attack  where  full  leading  stall  first  occurs.  A  further  increase  in 
angle  of  attack  results  in  a  decline  in  the  normal  force  coefficient. 
The  peak  value  of  the  normal  force  coefficient  is  greater  and  occurs  at 
a  higher  angle  of  attack  for  the  cascades  with  greater  stagger. 
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Abstrac  t 

Rotating  stall  In  a  compressor  stage  Is 
simulated  by  the  numerical  solution  of  a  fluid 
dynamic  model  of  an  Inlet  guide  vane  and  rotor 
row.  The  mathematical  formulation  of  the  model  Is 
described.  Solutions  for  the  flow  upstream  and 
downstream  of  the  modeled  stage  are  obtained  for 
normal  unstalled  operation  and  for  a  simulation  of 
stable  rotating  stall.  Examples  of  the  predicted 
velocity  and  pressure  fields  are  shown,  and 
conclusions  are  developed  regarding  recovery  from 
rotating  stall.  The  Influence  of  the  rotor  row 
loss  coefficient  curi>  and  other  blade  row  design 
variables  on  model  predictions  are  studied. 

Nomenclature 


BV  blade  velocity 

C  absolute  velocity 

F  Fourier  Transform 

F-1  Inverse  Fourier  Transform 

H  mechanical  energy  per  unit  mass  (fluid  head) 
i  blade  chord  length 

P  static  pressure 

P0  total  pressure 

q  volume  flow  rate 

s  blade  pitch 

t  time 

TF  transfer  function 

U  steady  axial  velocity 

u  axial  perturbation  velocity 

u  U  +  u 

V  steady  tangential  velocity 

v  tangential  perturbation  velocity 

v  V  +  v 

w  relative  velocity 


x  axial  coordinate 

y  tangential  coordinate 

C  vortlclty 

£  blade  stagger  angle  (ref.  axial  direction) 

p  density 

$  potential  function 

4  flow  coefficient 

w  total  pressure  loss  coefficient 

i  time  constant 

<1>  stream  function 

4  pressure  coefficient 

SUBSCRIPTS 

1  conditions  upstream  of  rotor 

2  conditions  downstream  of  rotor 

Introduction 

Axial-flow  compressor  post-stall  behavior  Is 
an  Important  aspect  of  aircraft  engine  stability 
and  operability.  The  rotating  stall  phenomenon, 
which  can  occur  under  certain  conditions,  limits 
the  range  of  engine  operation.  Reduced  compressor 
pressure  rise  and  mass  flow  associated  with 
rotating  stall  cause  performance  loss  and  engine 
overtemperature,  and  locally  unsteady  flow  within 
the  compressor  blading  can  produce  damage  through 
aeromechantcal  effects.  Compressor  recovery  from 
the  globally-stable  rotating  stall  condition  Is  an 
Important  concern. 

Conditions  which  give  rise  to  rotating  stall 
and  aspects  of  compressor  behavior  In  rotating 
stall  have  been  discussed  by  several  authors1'  '  . 
Analyses  by  Moore  have  shown  the  effect  of  the 
general  shape  of  the  compressor  characteristic  on 
post-stall  and  recovery  behavior.  The  present 
work  provides  predictions  of  compressor  stage 
post-stall  behavior  from  a  flow  model  which 
Includes  blade  row  design  features  and 
aerodynamics.  The  model,  based  on  an  ^pproach 
originally  developed  by  Takata  and  Nagano5,  shows 
local  features  of  rotating  stnll  In  a  compressor 
stage.  Including  flow  velocities  and  pressures. 
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Stall  Inception  and  recovery  points,  and  a  post¬ 
stall  operating  characteristic  are  predicted.  The 
details  of  the  loss  behavior  of  the  rotor  blade 
row  are  found  to  have  the  principal  Influence  on 
■odel  predictions. 


The  Model 

The  numerical  model  represents  components  of 
an  axial-flow  compressor,  that  Is,  a  set  of  Inlet 
guide  vanes  and  a  rotor.  In  the  model,  the 
compressor  rotor  Is  represented  by  a  discontinuity 
surface  (a  semi-actuator  disk)  at  which  the  energy 
transfer  due  to  the  rotor  blades  Is  Introduced 
lnco  the  flow.  The  flow  fields  upstream  and 
downstream  of  the  rotor  are  Interactively  coupled 
using  the  semi-actuator  disk.  The  model  Is  used 
to  predict  whether  the  onset  of  rotating  stall 
will  occur  under  any  particular  set  of  operating 
conditions  by  Imposing  a  perturbation  upon  the 
uniform  flow  field  of  the  upstream  face  of  the 
rotor  semi-actuator  disk.  If  the  rotor  Is 
operating  In  a  stable  regime,  the  perturbation 
velocity  will  diminish  with  time  and  the  flow  will 
return  to  Its  former  uniform  condition.  If  the 
rotor  flow  Is  aiade  unstable,  the  disturbance  will 
grow  Into  a  fully-developed  pattern  which  will 
propagate  around  the  rotor  seal-actuator  disk  In 
the  form  of  a  rotating  stall. 

The  solution  of  the  analytical  model  Is 
carried  out  using  numerical  techniques  to  solve 
the  unsteady,  nonlinear  equations  which  represent 
the  flow  fields  and  the  characteristics  of  the 
semi-actuator  disk. 

The  flow  fields  upstream  and  downstream  of 
the  rotor  are  represented  by  two-dimensional  flow 
models.  These  flow  fields  are  obtained  by 
unrolling  a  unit  thickness  cylindrical  section  of 
the  compressor  annulus  Into  a  two-dimensional 
plane,  as  shown  In  Fig.  1.  The  upstream  flow  Is 
assumed  to  be  Incompressible,  lnvlscld,  and 
lrrotatlonal.  The  downstream  flow  Is  considered 
Incompressible,  lnvlscld,  and  rotational. 
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In  models  of  this  type,  one  must  he  able  to 
calculate  the  dynamic  manner  In  which  the  blade 
row  characteristics  respond  to  changes  In  the 
inlet  flow.  The  present  stall  model  improves  on 
the  first-ordec  techniques  many  Investigators  have 
used  In  the  past  by  Incorporating  an  Improved 
dynamic  loss  model  developed  by  Sexton.  The 
model  utilizes  experimentally-measured  blade 
forces  to  describe  the  dynamic  blade  row  charac¬ 
teristics.  The  dynamic  loss  for  each  time  step  In 
the  numerical  calculation  Is  determined  by  first 
calculating  the  quasi-steady  variation  In  total 


pressure  loss  along  the  semi-actuator  disk 
representing  the  compressor  rotor,  employing  the 
quast-steady  loss  curve  and  the  flow  exit  angle 
for  the  rotor  cascade.  A  dynamic  response 
expression  in  the  form  of  a  transfer  function  (TF) 
Is  then  used  to  compute  the  dynamic  loss. 

The  Rotor  Semi-Actuator  Disk 


The  flow  upstream  of  the  rotor  semi-actuator 
disk  Is  assumed  to  consist  of  a  uniform  velocity 
field  on  which  a  perturbation  velocity  can  be 
superimposed.  Neglecting  body  forces,  the 
momentum  equations  for  the  perturbation  flow  field 
take  the  form  of  the  Fuler  Equations, 
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The  total  fluid  mechanical  energy  per  unit  mass  Is 

H  i(u2  +  v2)  •  (3> 

Taking  the  partial  derivative  with  respect  to  the 
x  and  y  directions  yields  for  the  perturbation 
flow  field 
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Combining  Eqs.  (1)  and  (4) 
lrrotatlonal  flow  field 
equations. 


(5) 


,  and  (2)  and  (5)  for  an 
yields  the  following 
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which  describe  the 

rotor. 

flow  field 

upstream  of 

the 

The 

upstream 

boundary  conditions 
perturbation  flow  field 

Imposed  on 
are  that: 

the 

(1) 

the  perturbation  velocity  vanishes  far  upstream, 
(2)  the  side  boundaries,  which  represent  a  common 
axial  trace  In  the  compressor  flow,  are  taken  as 
periodic  (l.e.,  any  flow  leaving  one  side  boundary 
must  enter  the  other  side  boundary  with  the  same 
velocity  vector),  and  (3)  the  boundary  condition 
at  the  rotor  face  is  a  prescribed  potential 
distribution  $(y)  at  the  Initial  time  which 
represents  the  perturbation  velocity.  A  new 
potential  distribution  4dy)  at  this  boundary  Is 
calculated  at  each  time  step.  The  solution  to 
Laplace’s  equation  with  this  boundary  condition 
yields  an  expression  for  the  axial  perturbation 
velocity  at  the  rotor  upstream  face. 


The  flow  field  downstream  of  the  rotor  can  be 
modeled  In  a  similar  manner  with  the  exception 
that  the  flow  can  no  longer  be  considered 
lrrotatlonal.  Vortlcity  Introduced  at  the  rotor 
semi-actuator  disk  will  he  convected  downstream, 


aaking  the  flow  rotational.  The  downstream  flow 
can  then  be  described  by  a  stream  function 
Kx.y.t)  and  a  vortlclty  distribution 
C(x,y,t).  The  vortlclty  Is  distributed  In  the 
flow  field  downstream  of  the  blade  row  and  the 
equations  of  aotlon  take  the  form  of 
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Combining  equations  (8),  (9),  (15),  and  (16) 

(9)  yields 


Elimination  of  H  from  equations  (8)  and  (9)  yields 
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+  (U  +  u)  -^  +  (V  +  v)  -  0  .  (10) 

In  developing  an  equation  to  represent  the 
exchange  of  energy  between  the  rotor  and  the  flow 
field,  the  time-dependent  energy  equation  for  a 
control  volume  Is  considered  In  the  form 

PQ(Hj  “  -  HiogJ)  ”  pa*  cosC(— 

(11) 
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u.  -  F_I  ITF  x  P(u)  ,  .  1 

dyn  quasi-steady 


(17) 


where  the  H^oss  term  Is  the  dynamic  value  of  the 
rotor  row  losses.  The  control  volume  considered 
for  Eq.  (11)  has  sides  parallel  to  the  flow  direc¬ 
tion  and  to  the  blade  pitch,  and.  Is  equal  to  the 
blade  pitch  In  width  along  the  circumferential 
direction.  The  terms  on  the  left  hand  side  of  the 
equation  represent  the  fluid  head  of  the  flow 
entering  the  control  volume,  the  fluid  head  of  the 
flow  leaving  the  control  volume,  and  the  total 
fluid  head  loss  In  the  control  volume  due  to 
viscous  effects.  The  terms  on  the  right  hand  side 
of  the  equation  represents  the  change  In  energy  In 
the  control  volume  and  the  change  In  position  of 
the  control  volume  as  It  moves  relative  to  the 

stationary  frame  of  reference.  Since  the  flow  Is 
assumed  Incompressible,  only  the  change  In  kinetic 
energy  of  the  flow  In  the  control  volume  need  be 
considered.  If  the  average  velocity  of  the  flow 
In  the  control  volume  Is  assumed  to  be  In  the 

direction  parallel  to  the  blade,  the  flow  will 

then  be  at  an  angle  £  to  the  axis  of  the 

compressor  and  the  time  variation  of  the  kinetic 
energy  In  the  control  volume  Is 


as  previously  discussed.  Equation  (17)  then  gives 
the  time-variation  of  Vj  for  the  rotor. 

The  Downstream  Flow  Field  Model 


The  two-dimensional  flow  field  downstream  of 
the  rotor  row  Is  assumed  to  be  Incompressible, 
lnvtscld,  and  rotational,  since  the  vortlclty  at 
the  seal-actuator  disk  representing  the  rocor  Is 
convected  downstream.  The  flow  downstream  of  the 
rotor  row  can  be  described  by  a  stream  func¬ 
tion  Kx,y,t)  and  a  vortlclty  distribution 
C(x,y,t).  The  equations  describing  the  down¬ 
stream  flow  field  are 


and 
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where 


u 


U  +  u 


The  position  variation  of  the  kinetic  energy  In 
the  control  volume  Is 


OHc.v. 


BV  £  (i  C2  sec2U 


sec2  7,  BV 


(11) 


Including  this  result  In  the  time-dependent  energy 
equation  yields: 
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v  *  V  ♦  v 

,  dv  _  du 
dx  dy 


The  solution  of  these  equations  Is  performed  l>y 
using  a  finite  difference  technique  throughout  the 
flow  field.  Since  u  and  v  are  known  leaving  the 
rotor  row,  the  vortlclty  can  be  calculated.  Far 
downstream  from  the  rotor  row,  the  boundary 
condition  Imposed  Is 


OP  _ 

dy 


0 


Taking  the  derivative  of  this  equation  In  the 
tangential  direction  yields  for  the  seml-actua tor 


(20) 


The  nuaerlcal  solution  of  the  downstream  flow 
field  is  performed  through  the  use  of  a  successive 
over  retaliation  technique. 


Once 

the  flow  field  Is 

determined,  a 

pressure 

solution 

If  performed  by 

solving 

the 

Poisson 

pressure 

equation 

.  2 

2  .2 
dx  by 

• 

(21) 

Resul ts 

The 

aodel  may  be 

util  lzed 

to 

develop 

solutions 

for  the  unsteady 

flow  fields 

ups  trean 

and  downstream  of  the  rotor 

actua tor 

disk 

•  Input 

data  are  the  quasi-steady  loss  and  turning  of  the 
modeled  rotor  blade  cascade,  as  determined  from 
cascade  tests  or  nusierlcal  calculations.  Figure  2 
shows  a  loas  curve  developed  from  available  data 
for  a  65  aeries  cascade  at  25°  stagger  angle.  In 
addition,  rotor  blade  velocity  BV,  blade  row 
dimensions  and  Inlet  air  density  are  specified. 


ANGLE  OF  ATTACK 


Fig.  2  Lots  Curve  for  W  Sorias  Coacode 

The  Initial  flow  field  for  the  model  is 
established  by  specifying  an  angle  of  attack  for 
the  flow  over  the  rotor  blades.  This  value,  with 
the  specified  geometric  parameters,  generates  a 
mean  axial  velocity  which  Is  maintained  throughout 
a  simulation  run.  The  solution  proceeds  through 
small  steps  In  time,  developing  perturbation 
values  of  velocities  (u  and  v).  If  an  Initial 
small  circumferential  velocity  perturbation  grows 
for  the  specified  angle  of  attack,  a  rotating 
stall  Is  simulated.  If  the  perturbation  becomes 
smaller,  the  specified  angle  of  attack  will  not 
produce  rotating  stall  In  the  model,  and  Is 
therefore  considered  to  he  less  than  the 
"Inception"  value. 

Since  the  model  calculates  local  velocity  and 
pressure  for  the  flow  field  In  the  region  of  the 
rotor,  the  variation  of  these  values  can  be 
examined  at  each  axial  station  In  the  flow  field. 


Sample  results  were  developed  for  a  blade  row 
having  loss  characteristics  as  shown  In  Fig.  2, 
and  for  an  initial  angle  of  attack  of  29*.  This 
was  well  Into  the  rotating  stall  region  for  this 
blade  row,  as  the  "inception"  point  was  found  In 
successive  runs  to  be  near  26°.  Figure  3  shows 
the  variation  in  axial  velocity  Immediately 
upstream  of  the  rotor  as  the  stall  develops  and 
become  stable.  Forty-nine  nodes  In  the 

circumferential  direction  were  provided  In  the 
nuaerlcal  formulation,  representing  0-2*  radians 
around  an  annulus.  The  Initial  imposed  perturba¬ 
tion  Is  shown  as  a  dashed  line,  while  the  final 
stable  rotating  velocity  field  Is  shown  as  a  bold 
solid  line.  Other  velocity  plots  in  Fig.  3 
represent  the  growing  perturbation,  and  the 
continuing  stable  axial  velocity  disturbance 
representing  a  rotating  stall.  The  simulated 
rotating  stall  Is  moving  at  a  rate  corresponding 
to  approximately  45X  of  the  blade  velocity. 
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The  local  angle  of  attack  of  the  flow  at  the 
rotor  face,  shown  in  Fig.  4,  Is  the  determining 
variable  for  the  look-up  of  quasi-steady  loss  and 
turning  angle  values.  The  dynamic  loss  term 

iVj  of  Eq.  (17)  Is  formed  through  the  use  of  a 
response  function. 

The  associated  tangential  or  circumferential 
velocity  field  immediately  upstream  of  the  rotor 

Is  shown  In  Fig.  5.  The  blockage  of  the  stalled 
rotor  passages  Induces  a  variation  In  tangential 

velocities  as  shown.  Figure  6  shows  the 

tangential  velocity  perturbations  downstream  of 
the  rotor  actuator  disk,  now  relatively  smaller 
because  of  the  swirl  velocity  Induced  by  the 
rotor.  The  varying  pressure  distribution  around 
the  simulated  annulus  downstream  of  the  rotor  Is 
shown  In  Fig.  7,  as  calculated  from  Eq.  (21).  The 
static  pressure  pattern  behind  the  rotor  arises  as 
a  consequence  of  the  varying  vortlctty  field 

behind  the  rotor.  The  downstream  pressure  distri¬ 
bution  changes  from  the  pattern  shown  In  Fig.  7  to 
a  uniform  pressure  far  downstream  of  the  rotor.  In 
accordance  with  the  condition  dP/dy  •  0,  Imposed 
by  Eq.  (20). 
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Employing  the  simulation.  It  Is  possible  to 
plot  axial  and  tangential  velocity  distributions, 
as  well  as  static  pressure  variations  at  any 
desired  point  in  the  modeled  flow  field,  upstream 
or  downstream  of  the  rotor.  The  plots  are  useful 
for  study  of  the  rotating  stall  phenomenon,  as 
well  as  for  comparisons  with  measured  values  from 
modeled  compressor  stages. 

Using  the  flow  field  pressures  and  velocities 
provided  by  the  model,  it  Is  possible  to  construct 
a  predicted  compressor  stage  characteristic  from 
the  results  of  the  simulation.  The  generally 
expected  form  of  the  rotating  stall  compressor 
characteristic  Is  shown  In  Fig.  8.  Study  of  the 
figure  shows  that  stable  operation  on  a  throttle 
line  within  the  bounded  region  Is  possible  at  two 
points.  The  actual  point  of  operation  depends 
upon  the  Immediate  past  history  of  the  compressor, 
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assuming  that  strong  externally-imposed  transients 
do  not  occur.  Thus,  a  compressor  In  stable  opera¬ 
tion  in  a  rotating  stall  will  remain  In  rotating 
stall,  so  long  as  the  throttle-imposed  flow 
coefficient  Is  within  the  range  of  the  rotating 
stall  characteristic.  Similarly,  a  compressor  In 
"normal"  operation  In  the  dual  stability  region 
will  remain  in  normal  operation  unless  the  flow 
coefficient  Is  reduced  below  the  "Inception 
value,  or  the  compressor  operation  is  upset  In 
some  way,  such  as  by  a  strong  transient  or  an 
Inlet  distortion.  The  existence  of  a  dual 
stability  region,  bounded  by  the  stall  Inception 
and  recovery  lines,  Is  the  essential  feature  of 
recovery  hysteresis.  8  measure  o(  recovery 
hysteresis  might  be  defined  as 


with  a  low  value  ot  k  representing  good  recovery 
characteristics. 


1 

I 


The  present  calculation  produces  an  operating 
characteristic  for  the  modeled  stage  which  Is 
similar  In  form  to  that  discussed  above.  It  will 
be  recalled  that  the  boundary  condition  far  down¬ 
stream  of  the  rotor  requires  a  uniform  static 
pressure  In  the  circumferential  direction.  Since 
the  local  velocities  at  this  point  are  also  known, 
a  mass-averaged  total  pressure  can  be  conveniently 
calculated.  The  total  pressure  rise  of  the  stage 
at  each  calculated  operating  point  can  thus  be 
expressed  In  the  form  of  a  pressure  coefficient, 

*-  ,  m> 

■j  p<Bvr 

where  the  station  numbers  refer  to  Fig.  1 


Specification  of  an  Initial  angle  of  attack. 
Inlet  guide  vane  discharge  angle,  and  rotor  blade 
velocity  produce  a  calculated  steady  axial 
velocity  U,  which  Is  maintained  as  a  convergence 
requirement  throughout  a  given  calculation.  Axial 
velocity  perturbations  (u)  developed  by  the  cosi- 
putatlons  are  added  to  the  constant  mean  axial 
velocity  (U)  to  obtain  the  local  absolute  axial 
velocity, 

u  -  u  +  U  .  (24) 

Defining  an  average  flow  coefficient 
49 
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Fig.  8  General  Form  of  Compressor  Characteristic 
•nth  Rotating  Stall. 


It  can  be  seen  that  a  jjtven  calculation  will  take 
place  at  a  constant  4.  Thus,  the  calculation 
simulates  the  development  of  a  rotating  stall  cell 
at  constant  mass  flow  or  f>,  rather  than  along  a 
fixed  throttle  line,  as  would  occur  In  a 
conventional  compressor  test.  A  series  of 
calculations  are  necessary  to  develop  the  rotating 
stall  operating  line,  each  defined  by  an  Initial 
angle  of  attack  of  the  rotor  flow.  The 
calculation  proceeds  from  the  uniform  axisymmetrlc 
flow  condition  at  Inception  along  a  line  of 
constant  4,  to  the  fully-developed  rotating  stall 
simulation. 

Figure  9  shows  a  stage  operating  characteris¬ 
tic  developed  from  a  series  of  eight  runs  of  the 
computer  model.  Again,  the  loss  curve  of  Fig.  2 
was  used,  and  a  blade  velocity  of  166  feet  per 
second  was  specified.  For  Initial  angles  of 
attack  of  25°  or  less,  the  Imposed  axial  velocity 
perturbation  decreased  In  magnitude,  predicting 
stable,  unstalled  operation.  For  an  Initial  angle 
of  attack  of  26”,  the  Initial  perturbation  grew  to 
a  final  stable  velocity  disturbance,  moving  at 
approximately  451  of  blade  velocity.  Following 
the  procedures  described  above,  the  values 
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of  ^  and  t  were  calculated  for  the  starting  point 
and  the  final  stable  rotating  stall  operating 
point.  Additional  calculations  at  29°,  31°,  and 
33°  Initial  angle  of  attack  were  performed, 
defining  the  predicted  post-stall  operating  line 
shown  in  Fig.  9.  The  starting  points  for  each 
calculation  define  a  line  of  axlsymmetrlc,  but 
unstable,  operation  of  the  stage  In  the  stalling 
regions  of  ♦. 

Constant  throttle  position  lines  are  sketched 
In  Fig.  9.  The  predicted  Inception  and  recovery 
throttle  positions  are  labeled,  leading  to  an 
estimate  of  the  recovery  parameter  R  -  0.05. 

Relatively  small  recovery  hysteresis  Is  evident 
here,  as  would  be  expected  for  the  simple  stage 
modeled.  As  can  be  scon  from  the  geometry  of  Fig. 
9,  the  use  of  Input  loss  curves  having  a  greater 
difference  between  the  minimum  and  maximum  loss 
values  would  lead  to  an  operating  line  In  stable 
rotating  stall  which  is  positioned  lower  on  a 
graph  such  as  Fig.  9.  Predicted  recovery 
hysteresis  In  the  form  of  the  recovery  parameter  R 
would  therefore  Increase. 

In  several  runs  of  the  model,  the  shape  of 
the  loss  curve  In  the  transition  region  from  low 
to  high  loss  has  been  shown  to  affect  the 

predicted  Inception  point  for  rotating  stall. 
Thus,  a  loss  curve  of  the  form  shown  In  Fig.  10 
would  lead  to  a  prediction  of  decreased  angle  of 
attack  for  stall  Inception,  a  less  abrupt  change 
In  the  operating  line,  and  lower  recovery 

hysteresis.  Figure  11  shows  a  stage  characteris¬ 
tic  generated  from  the  loss  curve  of  Fig.  10,  with 
all  other  parameters  held  constant.  Reduced 
vertical  spacing  of  the  pre-stall  and  post-stall 
operating  lines  Is  evident,  and  Improved  recovery 
behavior  Is  predicted.  Here,  the  R-parameter  Is 
calculated  as  essentially  zero,  quantifying  the 
Improved  recovery  prediction. 

The  model  Is  two-dimensional,  and  Is  limited 
In  prediction  accuracy  by  this  assumption.  The 
well-known  tip-stalling  behavior  of  a  single 


stage,  leading  to  a  "progressive",  rather  than  an 
abrupt  transition  to  rotating  stall,  is  not 
simulated.  The  predicted  behavior  Is  more  closely 
related  to  the  full  span  stalling  behavior 
observed  In  the  stages  of  multistage 
compressors.  Uhtle  three-dimensional  effects 
cannot  be  Included  In  the  present  formulation, 
extension  to  Include  local  reverse-flow 
calculation  capability  as  well  as  the  addition  of 
a  downstream  stator  row  Is  possible. 
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Fig.  10  Modified  Loss  Curve 

Conclusions 


A  mathematical  model  for  the  two-dimensional 
unstalled  and  stalled  flow  In  the  neighborhood  of 
a  compressor  stage  composed  of  an  Inlet  guide  vane 
row  and  a  rotor  row  has  been  described.  Principal 
input  data  for  the  model  are  the  rotor  row  quasl- 
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steady  loss  and  turning  angle  data,  a  dynamic  loss 
response  function,  the  blade  arrangement  In  terms 
of  stagger  and  chord,  and  the  rotor  blade 
tangential  velocity.  Predictions  of  local  velo¬ 
city  and  pressure  fields  In  rotating  stall  are 
obtained.  A  compressor  stage  performance  map  Is 
predicted.  Including  stable  pre-stall  and  post¬ 
stall  operating  lines,  and  a  measure  of  recovery 
hysteresis  is  suggested.  The  principal  Influence 
on  model  predictions  Is  found  to  be  contained  In 
the  shape  of  the  rotor  loss  curve. 
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